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ABSTRACT 

We use the ultra-deep WFC3/IR data over the HUDF and the Early Release Science WFC3/IR data 
over the CDF-South GOODS field to quantify the broadband spectral properties of candidate star- 
forming galaxies at z ~ 7. We determine the [/U-continuum slope /3 in these galaxies, and compare 
the slopes with galaxies at later times to measure the evolution in j3. For luminous L* =3 galaxies, we 
measure a mean t/V-continuum slope j3 of —2.0 ± 0.2, which is comparable to the (3 ~ — 2 derived at 
similar luminosities at z ~ 5 — 6. However, for the lower luminosity 0.1L* =3 galaxies, we measure a 
mean j3 of —3.0 ±0.2. This is substantially bluer than is found for similar luminosity galaxies at z^4, 
just 800 Myr later, and even at z~5-6. In principle, the observed of —3.0 can be matched by a very 
young, dust-free stellar population, but when nebular emission is included the expected /3 becomes 
> —2.7. To produce these very blue /3's (i.e., j3 ~ —3), extremely low metallicities and mechanisms 
to reduce the red nebular emission seem to be required. For example, a large escape fraction (i.e., 
fesc ^ 0.3) could minimize the contribution from this red nebular emission. If this is correct and 
the escape fraction in faint z ~ 7 galaxies is >0.3, it may help to explain how galaxies reionize the 
universe. 

Subject headings: galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

The spectral properties of high-redshift galaxies must 
undergo dramatic changes at some point in the past, 
as the metallicities in these systems drop to lower val- 
ues and these systems become progressively younger. 
In the limit of low metallicities, gas is no longer able 
to cool efficiently, likely resulting in massive extremely 
low-metallicity (or Population III) stars whose hot at- 
mospheres are expected to result in a very hard UV- 
continuum spectrum and strong Hell 1640 emission. 
However, the strong UV flux coming from hot stars is 
expected to be largely offset by the redder nebular con- 
tinuum light produced by the ionized gas surrounding 
these massive stars (e.g., Schaerer 2002). 

The newly installed WFC3/IR camera on the Hubble 
Space Telescope permits us to observe faint z > 7 galax- 
ies > 40 x more efficiently than before, providing us with 
our most detailed look yet at the UV light and spectral 
properties of z > 7 galaxies. Already >25 likely z^7- 
8 galaxies have been identified in the early WFC3/IR 
data over the Hubble Ultra Deep Field (HUDF: Oesch 
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et al. 2009b; Bouwens et al. 2009b; McLure et al. 2009; 
Bunker et al. 2009), and >20 z ~ 7 candidate galaxies in 
the WFC3 /IR Early Release Science (ERS) observations 
over the CDF-South (R.J. Bouwens et al. 2009, in prep). 

Here we take advantage of these early WFC3/IR ob- 
servations to study the spectral properties of candidate 
z ~ 7 galaxies. Our principal focus will be on the slope 
of z ~ 7 galaxy spectra in the Z7U-continuum - since 
this slope is the primary observable we can derive from 
the available broadband imaging data with WFC3. The 
C/U-continuum slope (3 (fx oc A' 3 : e.g., Meurer et al. 
1999) provides us with a powerful constraint on the age, 
metallicity, and dust content of high-redshift galaxies; it 
has already been the subject of much study at z ~ 3 — 6 
(Lehnert & Bremer 2003; Stanway et al. 2005; Yan et al. 
2005; Bouwens et al. 2006; Hathi et al. 2008; see Bouwens 
et al. 2009a for a systematic study at z ~ 2 — 6) and even 
at z ~ 7 (Gonzalez et al. 2009) using NICMOS data. 

Throughout this work, we quote results in terms of 
the luminosity L* z=3 Steidel et al. (1999) derived at z ~ 
3: M 17 ooab = -21.07. We refer to the HST F606W, 
F775W, F850LP, F105W, F125W, and F160W bands as 
Veoe, «775, z 850 , Yio5 5 ^125, and H 160 , respectively, for 
simplicity. Where necessary, we assume £Iq = 0.3, fl\ = 
0.7, Ho — 70km/s/Mpc. All magnitudes are in the AB 
system (Oke & Gunn 1983). 

2. OBSERVATIONAL DATA 

Our primary dataset is the ultra deep WFC3/IR data 
over the HUDF and the wide-area WFC3/IR ERS obser- 
vations over the CDF-South (PI O'Connell: GO 11359). 
The HUDF data permit us to identify lower luminosity 
z ~ 7 galaxies and quantify their properties (e.g., Oesch 
et al. 2009b), while the ERS data permit us to do the 
same for more luminous z ~ 7 galaxies. 

For our HUDF dropout selections, we make use of the 
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vl.O reductions of the HUDF ACS data (Beckwith et al. 
2006) rebinned on a 0.06"-pixel scale, and our own re- 
duction of the HUDF09 WFC3/IR data over the HUDF 
(Bouwens et al. 2009b; Oesch et al. 2009b). The op- 
tical ACS imaging over the HUDF reach to 29.4, 29.8, 
29.7, and 29.0 AB mag (5cr: 0.35"-diameter apertures) 
in the -B435, Veo6; 2775, #850 bands, respectively, while 
the near-IR WFC3/IR data reach to 28.8, 28.8, and 28.8 
(5c: 0.35" apertures) in the Y105, J125, and -H160 bands, 
respectively. The PSF FWHMs are ~0.10" for the ACS 
data and ~0.16" for the WFC3/IR data. 

For our dropout selections over the WFC3/IR ERS 
fields, we make use of our own reductions of the available 
ACS/WFC data over the GOODS fields (Bouwens et al. 
2006, 2007). Our reductions reach to 28.0, 28.2, 27.5, 
and 27.4 in the -B435, V e06 , z 775 , z 850 , respectively (5cr: 
0.35" apertures). This is similar to the GOODS v2.0 re- 
ductions, but reach ~0.1-0.3 mag deeper in the zsso-band 
due to the inclusion of the SNe follow-up data (Riess et 
al. 2007). Our reductions of the WFC3/IR ERS data 
was performed using the same procedures as we used on 
the HUDF09 WFC3/IR data. These data reach to 27.7 
and 27.4 in the J125 and Hieo bands, respectively (5a: 
0.35" apertures). 

3. RESULTS 

3.1. Catalog Creation 

Our procedure for source detection and photometry is 
similar to that used in previous studies by our team (e.g., 
Bouwens et al. 2007, 2009a) and relies upon the SExtrac- 
tor software (Bertin & Arnouts 1996) run in double image 
mode. Source detection is performed off the square root 
of x 2 image (Szalay et al. 1999: similar to a coadded im- 
age) constructed from all images redward of the Lyman 
break. Colors are measured in small scalable apertures 
with a Kron (1980) parameter of 1.2 (typically ~0.4"- 
diameter apertures). Fluxes measured in these small 
scalable apertures are then corrected (typically by ~0.4 
mag) to total magnitudes using the additional flux in a 
larger scalable aperture (Kron parameter of 2.5). 

3.2. z ~ 7 z-dropout selection 

For our lower-luminosity z ~ 7 z-dropout selection, 
we use a criterion very similar that used by Oesch et al. 
(2009b) over the HUDF09 WFC3/IR field: 

(^850 " ^105 > 0.8) A (Yi 05 - Ji 25 < 0.8) 

Over the ERS fields, no deep WFC3 Yios-band coverage 
is available, so we make use of the following criterion 
(R.J. Bouwens 2009, in prep): 

(2850 ~ J125 > 0.8) A ( J 125 - #i 6 o < 0.3) 

Sources are required to be undetected (< 2a) in all bands 
blueward of the break (and also undetected [< 1.5c] in no 
> 1 bands) to ensure that our selections are largely free 
of lower redshift contaminants. Sources must be detected 
at > 5. 5a in the J125 band to ensure that z ~ 7 candidates 
in our sample correspond to real sources. 

3.3. UV -continuum slope Measurements 

The C/U-continuum slopes /3 we estimate for sources in 
our z ~ 7 z-dropout sample are derived from the broad- 



-1 



-2 



-3 



-4 



-21 -20 -19 -18 

^UV.AB 

Fig. 1. — [/V-continuum slope j3 versus absolute UV magnitude. 
The red circles show the j3 determinations and absolute magnitudes 
(derived from the J125 - #160 colors and |(Ji25 + Hi6o)ab mag- 
nitudes, respectively) for individual z ~ 7 galaxies in our HUDF09 
and ERS selections. The large red squares (with lc error bars) 
represent the mean values in 1-mag bins {dotted lines). The black 
points correspond to the /3 determinations at z ~ 4 (Bouwens et 
al. 2009a). The /3's for lower luminosity galaxies at 2 ~ 7 (red 
circles) are much bluer (by A/3 of ~1) than those derived at z ~ 4 
(a 4c difference), (inset) The relative volumes available for select- 
ing galaxies with various /3's using our HUDF09 z-dropout crite- 
rion (see §3.4)). These volumes do not depend significantly on /3, 
demonstrating that the blue /3's observed for faint z ~ 7 galaxies 
is not a selection effect. 

band color J125 — #160 as 

/3 = 4.29( J125 ~ H 100 ) - 2.00 (1) 

The J125 and i?i6o bands here probe rest-frame 
~1550Aand ^1940A, respectively, for the typical z ~ 7 
z-dropout candidate in our sample, and are not affected 
by Lya emission or IGM absorption (in contrast to the 
Yi05-band). The above equation is derived assuming a 
flat spectrum source with no absorption lines (found to 
work very well by Meurer et al. 1999). 

In Figure [TJ we show the /3 determinations for z ~ 7 
z-dropout candidates versus luminosity. Both our ultra- 
deep HUDF09 and wide-area ERS selections are in- 
cluded, as are the z ~ 4 selections from Bouwens et al. 
(2009a). The trend of /3 with luminosity is illustrated 
with the large squares. The difference between z ~ 7 
and z ~ 4 is striking. 

Use of the photometry from Oesch et al. (2009b), 
McLure et al. (2009), or Bunker et al. (2009) for our 
z ~ 7 sample yield very consistent colors and similarly 
blue (or even bluer) /3's. In addition, use of a consistent 
0.7"-diameter aperture for the color measurements yield 
equally blue /3's (albeit with larger random uncertainties, 
due to the use of larger apertures). 

3.4. Consideration of Possible Selection Biases 

The distribution of /3's derived for sources in our z ~ 7 
sample can be affected by the selection process. This 
effect is seen in lower redshift samples, where galaxies 
with bluer /3's have larger selection volumes than galax- 
ies with redder /3's. To determine the importance of such 



LI' !..■■'••■' 



•• ■ z~4- ■• 



• — • • • 

z~7 



0.5 



> 0.3 I 

■3 2 
w 0.1 







-1 — 1 — 1 — I — r 



_L 



-4 -2 
J 1 1 1 1 l_ 



3 



-1 



1 1 1 1 1 1 1 1 1 
- Lu 

■ ^ i ^ 


1 i i i i 1 i i i i 1 i i i i 1 i i 
minous (L; =3 ) - 

! i i ^ 

x f 9- 

1 i i i i 1 i i i i 1 i i i i 1 i r 


i M M M M i M i 

— Lower Luminos 
- 1 i i i 1 i i i i 1 i i i i 


i i i i ii i i i ii i 

ity (o.i l;. 3 ) - 

1 1 • 

1 1 1 1 1 1 1 1 1 1 1 r 



2 3 4 5 6 7 
Redshift 

Fig. 2. — Mean tTV-continuum slope measured for galaxies at 
z ~ 7 (large red circles) versus galaxies of similar luminosities at 
lower redshifts. The top panel shows the evolution in this slope 
for galaxies with a UV luminosity of L*_„ (—21 AB mag) and the 
bottom panel shows this evolution for UV luminosities of 0.1L*_ 3 
(—18.5 AB mag). The determinations of Bouwens et al. (2009a) 
based upon ACS+NICMOS data are shown with the open red cir- 
cles. The short error bars are the random errors while the longer 
error bars also include possible systematic errors (e.g., see Bouwens 
et al. 2009a). We also show j3 determinations by Stanway et al. 
(2005: cyan square), Ouchi et al. (2004: blue pentagon), Adel- 
berger & Steidel (2000: black star), and Hathi et al. (2008: green 
open squares). Note the dramatic change in /3 in the ~800 Myr 
from z ~ 7 to z^4 for lower luminosity galaxies. 

effects, we constructed models with various /3 distribu- 
tions, added galaxies with these distributions to the ob- 
servational data, and then attempted to reselect these 
galaxies using the selection criteria in §3.2 to estimate 
the effective selection volume versus /3. We modelled the 
pixel-by-pixel profiles of the sources using similar lumi- 
nosity z ~ 4 £?-dropouts from the Bouwens et al. (2007) 
HUDF sample, but scaled in size (physical) as (1 + z)^ 1 
to match the observed size-redshift relationship (Oesch 
et al. 2009c; Ferguson et al. 2004; Bouwens et al. 2004). 

The selection volumes derived from these simulations 
are shown in the inset to Figure [T] Encouragingly, these 
volumes show only a modest dependence upon the input 
for the range —4 < (3 < — 1, and so selection biases are 
minimal. 

3.5. Low-redshift Comparison Samples 

To interpet the /3's we derive from our z ~ 7 z-dropout 
selections, we compare them with the /3's found for sim- 
ilar luminosity galaxies at z ~ 4 — 6. Bouwens et al. 
(2009a) provide determinations of these slopes as a func- 
tion of luminosity at z ~ 2 — 6. The /3 measurements at 
z ~ 4 are also shown in Figure []] for comparison. 



We can take advantage of the very deep, high-quality 
WFC3/IR data to obtain self-consistent determinations 
of these slopes at z~5-6. The z ~ 5 V and z ~ 6 i- 
dropouts over the HUDF09 and GOODS ERS data are 
selected in the same way as in Bouwens et al. (2007: 
see also Giavalisco et al. 2004; Beckwith et al. 2006), 
except that at z ~ 6 we also require galaxies to satisfy 
the criterion (zsso ~ J125 < 0.6). At z ~ 5 and z ~ 6, the 
/3's are estimated based upon the zgso — (^ios + Ji2s)/2 
and Y105 — (2J125 + i?i6o)/3 colors, respectively, which 
are a good match in rest-frame wavelength to the J125 — 
7?i6o colors used to estimate /3 at z ~ 7. The conversion 
formulas we use are 

P = 4.07(z 85 o - (Y W5 + Ji2 B )/2) - 2.00 (z ~ 5) (2) 
8 = 3.78(r 105 - (2 J 125 + JZibo)/3) - 2.00 (z ~ 6) (3) 

In Figure [2j we plot these /3 determinations as a func- 
tion of redshift, for both luminous L* galaxies and lower 
luminosity galaxies. 

3.6. z ~ 8 Comparison Sample 

Faint z ~ 8 Fios-dropout candidates in the HUDF09 
WFC3/IR data (Bouwens et al. 2009b) also have very 
blue J125 — -H16O colors, again suggesting /3's of ~ —3. We 
decided not to consider this sample here because the J125- 
band flux can be affected by Lyman series absorption and 
Lya emission at z > 8.1, and hence the results would be 
less robust. 

4. DISCUSSION 

Before discussing the extraordinarily blue UV- 
continuum slopes /3 ~ — 3 found for lower luminosity 
galaxies at z ~ 7, we first consider the relatively lumi- 
nous L* z=3 galaxy candidates. These sources have mea- 
sured /3's of —2, which is similar to that observed for 
luminous galaxies at z ~ 3 — 6 (Figure [2]). Such slopes 
can be fit by a moderately young, subsolar (0.2 Zq) stel- 
lar population, with a maximum E(B — V) of ~0.05 
(Calzetti et al. 2000 extinction law: see also Bouwens 
et al. 2009a). 

The lower luminosity (0.1L* =3 ) galaxy candidates at 
z ~ 7, by contrast, have observed /3's of —3.0 ± 0.2. 
This is much bluer than for luminous z ~ 7 galaxies (by 
A/3 ~ 1) and also much bluer than is found at z ~ 5 — 6 
(Figures [T][2]). This makes these galaxies of great interest 
since they are likely to be even younger, more metal poor, 
and dust-free than any galaxies known. 

This is not a selection effect (§3.4), and cannot be at- 
tributed to Lya emission contributing to the broadband 
fluxes. Lya does not move into the Ji25-band (used to 
estimate B) until z > 8.1, and < 10% of our z-dropout 
sample extends to z > 8. 

An AGN contribution seem similarly unlikely, given 
the rarety of AGN signatures in faint Lyman-Break and 
Lya-emitter samples (e.g., Nandra et al. 2002; Ouchi et 
al. 2008). 

What then is the explanation for the very blue /3's? 
We explore several possibilities: 

4.1. Standard stellar population models 

A first question is whether it is possible to obtain a 8 
of ~ —3 using standard stellar population models (e.g., 
Leitherer et al. 1999; Bruzual & Chariot 2003). The 
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Fig. 3. — i7V-continuum slope /3 we would expect for z ~ 7 galax- 
ies as a function of age for constant star formation (CSF) models 
and an instantaneous burst (Schaerer 2002). The gray band de- 
notes the observed mean /3 and its uncertainty. The slopes /3 de- 
rived from the stellar light (dotted lines) and the stellar + nebular 
light (solid lines) are shown. While it is in principle possible to 
obtain /3's of —3 with standard low mctallicity (>0.02 Zq) models, 
including the nebular emission associated with hot stars make the 
predicted /3's > —2.7 and hence too red. 
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Fig. 4. — {/ V-continuum slope /3 Schaerer (2003) calculated as a 
function of age for instantaneous burst models for different metal- 
licities. Both the slopes /3 derived from the stellar light (dotted 
lines) and the stellar + nebular light (solid lines) are shown. The 
pure stellar light (dotted lines) has very blue l/V-continuum slopes 
P's (/3 < —3) for all the low metallicity cases considered here. 
Changing the IMF does not appear to change the conclusion here 
in any significant way. Of course, these same hot low metallicity 
stars also ionize the gas around them, thus producing a substantial 
amount of redder nebular continuum light. This makes the total 
SED of a galaxy much redder in general, and in the calculations 
by Schaerer (2003) shown here, /3 never becomes bluer than —3.0. 



answer is that it is possible, but only for very young 
(<5 Myr) star-forming systems (see Figure [3]). However, 
to do so ignores the nebular continuum emission from 
the ionized gas around the young stars. Including this 
nebular continuum emission can redden the observed /3 
by as much as A/3 ~0.5. 

Figure [3] also shows the /3's predicted for several low 
metallicity (0.02 Zq and 0.2 Zq) starbursts, as a function 
of age for the Schaerer (2002) stellar population models 
(which - like Starburst99 Leitherer et al. 1999 - include 
a nebular contribution). In the best cases, the models 
predict /3's as blue as —2.7, which is redder than what 
we observe in our faint samples This would suggest 
that lower metallicities are needed since the standard 
Leitherer et al. (1999) or Bruzual & Chariot (2003) stellar 
population models do not produce blue enough colors to 
match those found in our lower luminosity z ~ 7 sample. 
Of course, the significance of this result is only modest 
(< 1.3er), but the similarly blue colors observed for z^8 
selections (Section 3.6) suggest that we may want to take 
this finding at face value. 

4.2. Extremely low metallicities (< 10~ 3 Zq) 

In Figure @] we present the /3's predicted by the 
Schaerer (2003) stellar population models (which con- 
veniently provide predictions at extremely low metallici- 
ties) as a function of age for instantaneous bursts assum- 
ing a metallicity of 0.5 x 10~ 3 Z Q , 0.5 x 10 _5 Z Q , and zero 
(population 111)0 From the figure, it is clear that the 
nebular component contributes significantly to the total 
light output from < 10 -3 Zq stellar populations. 

While initially somewhat red due to the nebular contri- 
bution, the predicted /3's for these ultra-low metallicity 
models become much bluer at ages > 10 Myr, eventually 
reaching /3's of —3. Such metallicities and ages are not 
necessarily unreasonable for lower luminosity galaxies at 
z ~ 7, and therefore at least one possible explanation for 
the very blue /3's in our selections is that the metallicities 
for galaxies in our sample may be < 10 -3 Zq. 

The above explanation may explain some of the very 
blue galaxies in our selection, but given that /3 ~ — 3 only 
for a limited period (10-30 Myr after a burst), it seems 
unlikely to be the general explanation (unless updated 
models revise the theoretical SEDs). 

4.3. Top-heavy IMF 

One seemingly attractive explanation for the very blue 
/3's observed is through a top-heavy IMF, since galaxy 
stellar populations would be weighted towards massive, 
blue stars. The difficulty with this explanation is that 
these same massive stars are extraordinarily efficient at 
ionizing the gas around them - resulting in substantial 
nebular emission and leaving the galaxy with a net ft no 
bluer (and likely redder) than the young (<1 Myr) bursts 
shown in Figures [SJH] (see also discussion in Leitherer & 
Hcckman 1995). 

9 To match the wavelength baseline for /3 measured here, we 
added -0.1 to the /3's (1300-1800Abaseline) tabulated in the 
Schaerer (2002, 2003). 

10 The use of instantaneous burst models allows us to explore 
the most extreme cases. Other models (r, inverse t, CSF) produce 
comparable /3's. 
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Fig. 5. — Mean [/V-continuum slopes /3's predicted for high- 
redshift galaxy samples versus the escape fraction fesc- The predic- 
tions are made using the Schaerer (2003) stellar population models 
(including nebular emission) for metallicities of < 10 — 3 Zq. Galax- 
ies are assumed to be observed at some random point during their 
star formation histories and to experience instantaneous bursts of 
star formation every 10-40 Myr (the dashed, solid, and dotted lines 
give the average fl for the first 10, 20, and 40 Myr, respectively, of 
the instantaneous burst). For fesc = 1, all of the ionizing radiation 
from a galaxy escapes into the IGM and hence does not contribute 
to ionizing the gas within a galaxy (and hence the contribution 
from nebular continuum emission to the total light is minimal). 
On the other hand, for fesc ~ (preferred in some simulations: 
e.g., Gnedin et al. 2008), the ionizing radiation from the hot stars 
does not make it out of galaxies — resulting in substantial nebu- 
lar continuum emission (see Figure [1} and hence a much redder 
0. Perhaps the very blue /3's observed (shaded gray region) could 
indicate that the escape fraction is larger at fesc > 0.3 than what 
has commonly been considered? 

4.4. Minimizing nebular emission 

A significant obstacle to matching the very blue /3's 
observed is the red nebular emission associated with hot, 
ionizing stars. The nebular contribution could be re- 
duced in a number of ways, by changes to the ioniza- 
tion parameter, metallicity, geometry, etc. Assessing the 
impact of such changes would benefit from further de- 
tailed modelling. However, we should emphasize that 
regardless of changes to the nebular contribution very 
low metallicity models (or very young ages) appear to be 
needed to match the very blue /3's observed. 



4.5. Changes in Escape Fraction? 

One possibility that we have explored to reduce the 
nebular emission contribution is if the ionizing radiation 
leaks directly into the IGM (e.g., due to the effect of SNe 
on the galaxies' ISM, possibly from a top-heavy IMF: 
Trenti & Shull 2009). We consider such a possibility 
schematically in Figure [5j showing the time-averaged /3's 
expected for stellar populations of various metallicities as 
a function of the escape fraction f esc of ionizing photons 
into the IGM. For f esc of unity, we simply recover the /3's 
from the pure stellar SEDs and for f esc of zero, we recover 
the stellar + nebular SEDs. For fractional f esc , we in- 
terpolate between the two extremes. The time-averaged 
/3's used for Figure [5] is based upon those presented in 
Figure 01 

Comparing the predicted /3's with that observed (grey- 
shaded region: Figure [5]), we see that f esc >0.3 would 
permit us to easily match the observations. Such an es- 
cape fraction is signicantly higher than the ~10% fre- 
quently assumed in calculations assessing the sufficiency 
of galaxies to reionize the universe. Since most of the 
luminosity density at z > 7 comes from low luminosity 
galaxies, the estimated number of ionizing photons in the 
z > 7 universe could increase by factors of >3. Such a 
large change could provide the needed photons to reion- 
ize the universe, providing a resolution to the current 
debate (e.g., Bouwens et al. 2008; Oesch et al. 2009a, 
2009b; McLure et al. 2009; Bunker et al. 2009; Gonzalez 
et al. 2009; Ouchi et al. 2009; Pawlik et al. 2009). 

5. SUMMARY 

The strikingly blue UV-continum slopes /3's seen at 
lower luminosities in candidate z ~ 7 galaxies (also ap- 
parent in the Bouwens et al. 2009b z ~ 8 sample) indicate 
that we are now beginning to explore a regime where the 
nature of galaxies and their stellar populations are under- 
going a dramatic change. These results raise many as yet 
unanswered questions, but could be heralding the transi- 
tion from the first stars and youngest objects within the 
first 400 Myr at z >10 to the galaxies that dominate the 
universe for the next 2 Gyr and provide clues as to the 
source of photons that reionize the universe. 



We thank Daniel Schaerer for helpful conversations. 
We are grateful to all those at NASA, STScI and 
throughout the community who have worked so diligently 
to make Hubble the remarkable observatory that it is to- 
day. We acknowledge the support of NASA grant NAGS- 



REFERENCES 



Adelberger, K. L. & Steidel, C. C. 2000, ApJ, 544, 218 

Beckwith, S. V. W., et al. 2006, AJ, 132, 1729 

Bertin, E. and Arnouts, S. 1996, A&AS, 117, 39 

Bouwens, R. J., Illingworth, G. D., Blakeslee, J. P., Broadhurst, 

T. J., & Franx, M. 2004, ApJ, 611, LI 
Bouwens, R.J., Illingworth, G.D., Blakeslee, J. P., & Franx, M. 

2006, ApJ, 653, 53 
Bouwens, R. J., Illingworth, G. D., Franx, M., & Ford, H. 2007, 

ApJ, 670, 928 

Bouwens, R. J., Illingworth, G. D., Franx, M., & Ford, H. 2008, 
ApJ, 686, 230 

Bouwens, R. J., et al. 2009a, ApJ, 705, 936 

Bouwens, R.J., et al. 2009b, ApJ, in press, arXiv:0909.1803 
Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 



Bunker, A., et al. 2009, MNRAS, in press, arXiv:0909.2255 
Calzctti, D., Armus, L., Bohlin, R. C, Kinney, A. L., Koornneef, 

J., & Storchi-Bergmann, T. 2000, ApJ, 533, 682 
Ferguson, H. C. et al. 2004, ApJ, 600, L107 
Giavalisco, M., et al. 2004a, ApJ, 600, L93 

Gnedin, N. Y., Kravtsov, A. V., & Chen, H.-W. 2008, ApJ, 672, 
765 

Gonzalez, V., Labbe, I., Bouwens, R., Illingworth, G., Franx, M., 
Kriek, M., Brammer, G. 2009, ApJ, submitted, arXiv:0909.3517 
Hathi, N. P., Malhotra, S., & Rhoads, J. E. 2008, ApJ, 673, 686 
Kron, R. G. 1980, ApJS, 43, 305 
Lehnert, M. D. & Bremer, M. 2003, ApJ, 593, 630 
Leitherer, C, & Heckman, T. M. 1995, ApJS, 96, 9 
Leitherer, C, et al. 1999, ApJS, 123, 3 



6 



McLure, R., et al. 2009, MNRAS, in press, |arXiv:0909. 2437] 
Meurer, G. R., Heckman, T. M., & Calzetti, D. 1999, ApJ, 521, 64 
Nandra, K., Mushotzky, R. F., Arnaud, K., Steidel, C. C, 

Adelberger, K. L., Gardner, J. P., Teplitz, H. I., & Windhorst, 

R. A. 2002, ApJ, 576, 625 

Oesch, P. A., et al. 2009a, ApJ, 690, 1350 

Oesch, P.A., et al. 2009b, ApJ, in press, arXiv:0909.1806 
Oesch, P.A., et al. 2009c, ApJ, in press, arXiv:0909.5183 
Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713 
Ouchi, M., et al. 2004, ApJ, 611, 660 
Ouchi, M., et al. 2008, ApJS, 176, 301 
Ouchi, M., et al. 2009, ApJ, 706, 1136 

Pawlik, A. H., Schaye, J., & van Scherpenzeel, E. 2009, MNRAS, 
394, 1812 



Riess, A. G., et al. 2007, ApJ, 659, 98 
Schaerer, D. 2002, A&A, 382, 28 
Schaerer, D. 2003, A&A, 397, 527 

Stanway, E. R., McMahon, R. G., & Bunker, A. J. 2005, 

MNRAS, 359, 1184 
Steidel, C. C, Adelberger, K. L., Giavalisco, M., Dickinson, M. 

and Pettini, M. 1999, ApJ, 519, 1 
Szalay, A. S., Connolly, A. J., & Szokoly, G. P. 1999, AJ, 117, 6 8 
Trenti, M., & Shull, M. 2009, ApJ, submitted. larXiv: 0905. 4505 
Yan, H., et al. 2005, ApJ, 634, 109 



